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Abstract. We consider the effects of the Galactic magnetic 
field on the propagation of ultra high energy cosmic rays 
(UHECRs). By employing two methods of trajectory sim- 
ulation, we investigate the possibility that UHECRs are pro- 
duced within the Galaxy with paths strongly influenced by 
the Galactic magnetic field. Such trajectories have the poten- 
tial to reconcile the existing conflict between proposed local 
sources and isotropic UHECR arrival directions. 



1 Introduction 

As is well known, the Greisen-Zatsepin-Kuzmin (GZK) cut- 
off constrains detected UHECRs to have been produced in 
or near the Galaxy. Specifically, detected proton primaries 
with energies exceeding 5 x 10^^ eV must have been pro- 
duced within 50 Mpc of Earth, and nuclei propagation dis- 
tances are constrained even further (Puget et al., 1976). The 
near-isotropy of detected UHECR arrival directions, how- 
ever, suggests that Galactic source locations are not easily 
associated with the observed arrival directions. We attempt 
to reconcile these observations by examining the possibil- 
ity that UHECRs are Galactic in origin, but consist of iron 
nuclei primaries with trajectories influenced by the Galactic 
magnetic field. It has been shown by Blasi et al. (2000) that 
MHD winds from young neutron stars are capable of Galac- 
tic production of iron primaries that fit the observed UHECR 
energy spectrum. Such a production mechanism implies that 
most source locations should be found within the Galactic 
disk. With this in mind, we propagate iron nuclei UHECRs 
through a realistic model of the Galactic magnetic field to 
investigate the possibility that Galactic sources and isotropic 
arrival directions are not mutually exclusive phenomena. 
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2 Galactic Field Model 

Following Stanev (1997) and Harari et al. (1999), we adopt 
a large-scale regular Galactic magnetic field associated with 
the spiral arms of the Galaxy. Specifically, we choose a 
bisymmetric even-parity field model (BSS-S) in which the 
field reverses direction between different spiral arms, but is 
symmetric with respect to the Galactic plane. The field strength 
in the plane, directed along the spiral arms, at a point (p, 9) 
in Galactocentric coordinates is given by 

B3p = Bo(p)cos(0-/31n(p/po)) (D 

with po — 10.55 kpc as the Galactocentric distance to max- 
imum field strength at / = 0° with (3 = 1/ tanp, where the 
pitch angle is p = —10°. The radial dependence of the field 
strength is given by 

Bq{p) ^—t&n\i^{ — ) iiG (2) 
P Pi 

where tq = 8.5 kpc is the Galactocentric distance to the sun 
and pi = 2 kpc is a smoothing factor that allows for 1/p 
behavior beyond 4 kpc from the Galactic center The field 
equation in the Galactic plane is given by 

B{p,0,z = Q) = Bsp[sinpp + cospO] (3) 

For this even-parity model, we introduce a z dependence of 
the following form 

Bsip, 0, z) = Bip, 0,z^0) (^^J^ + ,(4) 

where the values zi = 0.3 kpc and 22 = 4 kpc reflect the 
scale heights of the field in the Galactic disk and halo, re- 
spectively. Figure (1) gives a graphical representation of this 
regular field model, clearly illustrating the field association 
with the spiral arms. Note that this model generates field val- 
ues ranging approximately from — 6 pG, with a value of 

3 pG in the solar neighborhood. 
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Fig. 1. Vector plot of the Galactic magnetic field. The solar position 
is shown at the intersection of the Galactic coordinate axes. 



3 Simulations 

We develop two distinct approaches to the simulation of Galac- 
tic UHECR trajectories. In our first method, we model a dis- 
tribution of Galactic sources and emit UHECRs according to 
an assumed N{E) oc emission energy spectrum in or- 
der to study the arrival energy spectrum at a detector. This 
allows for comparison between the injected and observed 
energy spectrum. Our second approach consists of model- 
ing the antiparticle trajectories as they depart Earth, given 
the detected UHECR arrival spectrum. This is equivalent to 
plotting particle trajectories that are guaranteed to intersect 
Earth, and this method has the advantage of allowing us to 
investigate potential source locations corresponding to real 
observed arrival directions. 

3 . 1 Particle traj ectories 

Our simulation of UHECR particle trajectories proceeds from 
a reasonable distribution of neutron stars within the Galactic 
disk to serve as injection regions for UHECRs. For simplic- 
ity, we assume that all sources are located in the Galactic disk 
of radius 25 kpc and thickness .65 kpc (taken from the scale 
height of the thin disk). We can derive a reasonable number 
of distinct sources from 

Nr^s = ^ (5) 

where tint = 30 Myr is the program integration time and 

fns — (1/100 years) is the neutron star birth rate. These 
values produce Nns = 300, 000 distinct UHECR sources. 
Given this set of sources, we assign a random time temit < 
tint at which each source emits randomly directed UHECRs 
according to the £^~^ energy spectrum, as is expected from 



neutron star sources (Blasi et al., 2000). The particles then 

propagate through the Galaxy with paths influenced by the 
Lorentz force as they traverse the Galactic magnetic field. 

Since an Earth-sized detector is too small to detect a sig- 
nificant number of events with computationally reasonable 
numbers of injections, we develop a series of larger detectors. 
The first is a 2D Galactocentric cylindrical detector of radius 
ro = 8.5 kpc, spanning 20-40 pc above the Galactic plane. 
The second detector is a 2D Galactocentric ring located 30 
pc above the Galactic plane with an inner radius of 8.49 kpc 
and an outer radius of 8.5 1 kpc. These 2D detectors, placed at 
the solar distance from the Galactic center, possess a number 
of advantages over local 3D detectors by generating a much 
higher detection flux while simultaneously reducing the bias 
for detection of local sources. 

With these detectors defined, we sample UHECR energy 
spectra to compare the detected energy spectrum with the as- 
signed injection energy spectrum. To accomplish this 
task, we inject millions of particles from our 300,000 sources 
and plot the detected energy spectrum. Figure (2) shows a 
typical detection energy spectrum from such a simulation. 
Since the energy emission spectrum is continuous, we bin 
the detected data before fitting it to a function of the form 
N{E) oc i?^". For the data shown, we calculate a = 
1.00 ± .09. Thus, we find that the emission and arrival spec- 
tra are not significantly different for UHECRs. This is an 
important result since it is unknown a priori if the detection 
energy spectrum should reflect the emission spectrum when 
particles have the potential to be trapped in the Galactic field. 



Fig. 2. Detected energy spectrum from particle simulations. Fit is 



of theformAf(S) oc 



, matching the emission spectrum. 



Fig. 3. 50 EeV antiparticle trajectories, emitted from Earth. The 
Galactic field is shown in the background for reference. 

3.2 Antiparticle trajectories 

This simulation addresses the other end of the UHECR prob- 
lem, propagating anti-iron nuclei from a distribution of ar- 
rival directions at Earth. In this approach, we assign energies 
to the emitted antiparticles and propagate them through the 
Galaxy to determine possible source locations for the corre- 
sponding UHECR particles. Figure (3) shows one such sam- 
ple set of trajectories for the case of isotropically distributed 
50 EeV iron nuclei. This method guarantees that our trajec- 
tories intersect Earth, but we are now left to evaluate source 
location without an initial neutron star distribution. Keep- 
ing in mind the need for Galactic neutron stars as sources, 
we define a potential UHECR source location as any point 
along an antiparticle trajectory that intersects the Galactic 
disk. Furthermore, we closely examine those antiparticle tra- 
jectories for which the path first leaves and then reenters the 
disk. These particular paths are the least likely to be easily 
identified with their Galactic sources, and their abundance 
would lend support to Galactic UHECR origins. 

We proceed by first emitting a known distribution of anti- 
iron nuclei from Earth. Specifically, we choose to trace back 
AGASA events with detected energies that exceed 100 EeV. 
Then, using the AGASA energies and arrival directions, we 
send the particles back through the Galactic field. Since most 
of the detected high-energy AGASA events are observed to 
have arrived away from the plane of the disk, we don't have 
to worry about the exact difference between arrival direction 
and source location. Any path that manages to reintersect the 
disk some time after emission will have deviated enough to 
be of interest to us. 

Figure (4) shows the set of trajectories modeled from those 



Fig. 4. Computed trajectories for AGASA events exceeding 100 
EeV with a 3 //G Galactic field. 

AGASA events above 100 EeV. It is clear that there is lit- 
tle turning back towards the Galactic disk once the particles 
have exited. The only path that remains in the Galactic disk 
for an extended period of time corresponds to the AGASA 
particle detected nearest the disk. With our choice of field 
parameters, the AGASA events do not clearly point back to 
sources within the Galactic disk. 



4 Discussion and conclusions 

Needless to say, the results of both types of simulation ex- 
hibit great sensitivity to certain parameters of the regular 
magnetic field model. Specifically, we find that variations 
in field strength strongly affect possible source locations for 
detected UHECRs. Figure (5) shows a contrast between the 
source location distributions for local field values of 3 and 
6 /xG. The 3 /zG field results in an abundance of detections 
involving local sources, while the 6 /iG case exhibits pro- 
nounced non-local Galactic features in the detection spec- 
trum. Furthermore, differences in the field strength do greatly 
alter our AGASA trajectories. Figure (6) shows the trajecto- 
ries produced from AGASA information with a 6 /iG Galac- 
tic field. Clearly, our choice of field parameters is important. 

In addition to the regular Galactic field used in our models, 
there are a number of additional magnetic field components 
that may be present in the Galaxy. In particular, a Galactic 
wind similar in nature to the solar wind can significantly alter 
UHECR trajectories if the wind field strength is comparable 
to the regular field strength. Such a wind model with an az- 
imuthal field strength of 7 fiG has been proposed to redirect 
the highest energy particles such that a single source in Virgo 
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Fig. 6. Computed trajectories for AGASA events exceeding 100 
EeV witli a 6 /nG Galactic field. 

can be detected as an isotropic distribution on Earth (Ahn et 
al., 1999; Harari et al., 2000). This strong wind effectively 
funnels isotropic arrival directions to the north Galactic pole 
for a narrow range of energies. For iron nuclei, this strong 
version of the Galactic wind guarantees diffusive behavior 
up to the highest energies detected thus far. 

Local high-field regions within the Galaxy also have the 
potential to contribute to the diffusion of UHECRs. Molec- 
ular clouds, for instance, are typically associated with mG 



fields that could serve as scattering regions for cosmic rays. 
Our preliminary research has indicated that the known distri- 
bution of molecular clouds has a small enough filUng factor 
to cause little change in the paths of most UHECRs, but these 
and other magnetic field inhomogeneities must be better un- 
derstood before the picture of UHECR propagation can be 
made complete. 
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